home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / slasq1.z / slasq1
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAASSSSQQQQ1111((((3333FFFF))))                                                          SSSSLLLLAAAASSSSQQQQ1111((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLASQ1 - SLASQ1 computes the singular values of a real N-by-N bidiagonal
  10.      matrix with diagonal D and off-diagonal E
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SLASQ1( N, D, E, WORK, INFO )
  14.  
  15.          INTEGER        INFO, N
  16.  
  17.          REAL           D( * ), E( * ), WORK( * )
  18.  
  19. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  20.         SLASQ1 computes the singular values of a real N-by-N bidiagonal
  21.         matrix with diagonal D and off-diagonal E. The singular values are
  22.         computed to high relative accuracy, barring over/underflow or
  23.         denormalization. The algorithm is described in
  24.  
  25.         "Accurate singular values and differential qd algorithms," by
  26.         K. V. Fernando and B. N. Parlett,
  27.         Numer. Math., Vol-67, No. 2, pp. 191-230,1994.
  28.  
  29.         See also
  30.         "Implementation of differential qd algorithms," by
  31.         K. V. Fernando and B. N. Parlett, Technical Report,
  32.         Department of Mathematics, University of California at Berkeley,
  33.         1994 (Under preparation).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      N       (input) INTEGER
  38.              The number of rows and columns in the matrix. N >= 0.
  39.  
  40.      D       (input/output) REAL array, dimension (N)
  41.              On entry, D contains the diagonal elements of the bidiagonal
  42.              matrix whose SVD is desired. On normal exit, D contains the
  43.              singular values in decreasing order.
  44.  
  45.      E       (input/output) REAL array, dimension (N)
  46.              On entry, elements E(1:N-1) contain the off-diagonal elements of
  47.              the bidiagonal matrix whose SVD is desired.  On exit, E is
  48.              overwritten.
  49.  
  50.      WORK    (workspace) REAL array, dimension (2*N)
  51.  
  52.      INFO    (output) INTEGER
  53.              = 0:  successful exit
  54.              < 0:  if INFO = -i, the i-th argument had an illegal value
  55.              > 0:  if INFO = i, the algorithm did not converge;  i specifies
  56.              how many superdiagonals did not converge.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.